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INTRODUCTION 


This  report  is  the  final  scientific  report  covering  the  period 
15  July  1987  thru  14  July  1990. 

In  the  last  research  Progress  and  Forecast  report  we  have 
developed  three  main  subjects: 

1  -  Propagation  of  regional  phases  in  western  Europe 

2  -  Far  field  effects  due  to  complex  geological  structures 

vicinity  of  point  sources. 

3  -  Source  evaluation  at  teleseismic  distance. 


1-  The  main  purposes  of  the  subject  n“  1  have  been  achieved  and 
we  present  a  complete  report  on  that  subject  including  our 
results  on  QLg  and  Q  coda  evaluations  for  a  set  of  regional 
events  recorded  on  an  homogeneous  seismic  network. We  found  very 
similar  results  in  the  evaluation  of  both  Q  factors  in  the 
studied  frequency  zone  (  l  to  15  Hz  )  consistent  with  other 
studies,  (e.g.:  Singh  and  Herrmann  1983)  and  suggesting  that 
both  QLg  and  Q  coda  measure  the  total  Q  factor  resulting  from 
intrinsic  and  scattering  attenuations. 

Concerning  the  interpretation  of  Q-values,  their  correlation 
with  the  tectonic  activity  has  been  established  in  France 
(Chapter  1).  As  a  complementary  study,  we  have  developed  a 
research  programm  on  the  attenuation  of  Lg  phases  propagating 
through  the  western  Alps  (Chapter  2).  We  have  collected  a  set 
of  records  from  short  period  seismic  stations  in  France  and 
Italy  for  earthquakes  which  occured  within  and  around  the  Alps. 
Lg  to  Pn  amplitudes  have  been  computed  for  different  frequency 
ranges  to  characterize  the  relative  efficiency  of  wave 
propagation  in  the  crust  and  in  the  upper  mantle.  Along  some 
paths  the  extinction  of  Lg  waves  is  almost  complete. 

From  this  computation,  we  have  built  up  a  map  showing  a  very 
sharp  zone  of  anomaly  of  propagation  located  on  the  positive 
Bouguer  anomaly  which  extends  along  a  narrow  band  in  the 
western  Alps. 

At  a  frequency  of  2Hz  the  value  of  the  mean  ratio  between  Lg 
and  Pn  amplitudes  is  larger  than  4  in  the  western  Alps  and  in 
central  France,  but  decreases  to  less  than  0.3  in  the  Bouguer 
anomaly  zone  or  to  the  east  of  this  zone. 

In  a  second  part,  we  have  modeled  the  propagation  of  Lg  waves 
through  the  Alps.  For  that  purpose  we  buixt  up  a  simplified 
cross  section  of  the  western  Alps  including  the  most  important 
features  likely  to  severely  affect  the  crustal  wave  propagation 
(geometry  f*oho;  deep  redi-ent  •''•'•y  basin  on  both  sides  of 

the  range ) . 
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The  computed  sismograms  obtenaid  by  the  discrete  wavenumber 
boundary  integral  methou  for  irregularly  layered  medium 
(Bouchon  1989)  limited  to  the  case  of  Sh  waves,  show  a  large 
decrease  of  Lg  amplitude  in  the  central  part  of  the  mountain 
range.  The  amplitudes  strongly  increase  when  the  waves  reach 
the  molassic  sedimentary  basin  to  the  west  of  the  Alps. 

Two  causes  might  account  for  the  results  :  the  curvature  of  the 
Moho  in  its  deepest  part  tending  to  focus  reflecting 
shear  waves  which  make  up  the  Lg  wavetrain:  the  diffraction  of 
the  crustal  waves  near  the  edge  of  the  molassic  basin  will 
result  in  the  excitation  of  waves  trapped  in  the  sediments 
associated  with  high  amplitudes  (Campillo  1987). 

2-  The  subject  n*  2  concerns  the  evaluation  of  farfield  effects 
due  to  complex  geological  structures  in  the  vicinity  of  point 
sources  ( Chapter  3 ) . 

Using  synthetic  models,  we  investigate  and  try  to  separate  the 
relative  effects  of  the  structure  heterogeneities  and  of  the 
source  medium  which  contribute  to  the  teleseismic  ground  motion 
recorded  by  a  thin  aperture  network.  The  numerical  simulation 
are  carried  out  using  the  discrete  wavenumber  -  boundary 
integral  equation  formulation.  The  computed  seismograms 
correspond  to  the  P-wave  displacement  component  recorded  at  a 
distance  of  80°.  In  the  calculation  the  source  region  is 
heterogeneous,  the  travel  in  the  mantle  is  assumed  to  be 
governed  by  the  Futterman  inelastic  attenuation  law,  and  the 
receiver  transfer  function  is  assumed  to  be  flat.  We  show  that 
the  shape  of  the  basement  has  a  large  influence  on  the 
magnitude  estimates.  The  magnitude  variations  related  to  the 
explosion  locations  can  get  as  large  as  mb  =  0.25,  and  the 
frequency  corresponding  to  the  maximum  observed  amplitude  is 
correlated  to  the  thickness  of  the  alluvium  deposit. 

This  study  shows  the  importance,  at  teleseismic  distances,  of 
the  structure  heterogeneities  surrounding  the  source  region, 
and  describes  the  evolution  of  the  seismogram  shapes  with  the 
complexity  of  the  source  model . 

3-  The  subject  n°3  deals  with  an  evaluation  of  seismic  sources 
at  teleseismic  distances  (chap.  4). 

The  main  results  were  already  presented  in  a  "research  progress 
and  forecast  report”  (11-02-1989)  but  we  think  worthwile,  for 
completeness,  to  summary  the  results  in  this  final  scientific 
report . 

The  records  obtained  by  the  French  Seismic  Network  on  about 
twenty  Eastern  Kazakh  events  wore  studied. 

The  first  goal  of  this  work  was  to  improve  the  magnitude 
determination.  Of  course,  the  relatively  small  aperture  of  the 
seismic  array  usad,  limits  this  determination.  The  only  hope 
was  to  improve  the  relative  magnitude  between  the  events,  as 
all  sources  are  in  the  same  small  geographic  area. 
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We  definitively  improve  the  magnitudes  of  those  events,  but 
only  slightly.  The  main  advantage  of  our  proposition  is  that 
this  determination  is  independant  of  the  seismologist  and  the 
main  disadvantage  is  that  it  has  to  be  calibrated,  so  is  only 
possible  to  apply  for  a  known  "seismic"  area. 

To  explain  the  observed  limitations,  we  look  at  the  time 
history  of  each  sources,  as  seen  from  the  French  Seismic 
Network.  The  already  known  division  of  the  Shagan  River  Test 
Site  in  a  SW  and  a  N£  sub  Test  sites  was  found  again. 

We  interprete  the  variation  in  the  seismic  sources  as  the 
crossing  by  seismic  rays  originated  in  the  N.E.  part  of  Shagar 
River  of  a  low  velocity  zone  in  the  vicinity  of  the  Test  Site. 

Certainly  we  are  close  from  the  limits  of  seismic  sources 
interpretations  using  a  small  aperture  seismic  array  at 
teleseismic  distances. 


Regional  seismic  phases  attenuation  and  seismic  source  evaluation 


J.L.  Plantet 

Laboratoire  de  Detection  et  de  Geophysique  Commissariat  Jt  l’Energie  Atomique  BP  12, 
91680  Bruy£res-le-Chatel,  France 

Y.  Cansi 

Laboratoire  de  Detection  et  de  Geophysique  Commissariat  il  1’Energie  Atomique  BP  12, 
9168C  Bruycies-I^-Chate',  France 

M.  Campillo 

Laboratoire  de  Gdophysique  Interne  et  Tectonophysique  Universite  Joseph  Fourier  and 
Observatoire  de  Grenoble  IRIGM,  BP  53X,  38041  Grenoble,  France 

1-7.1  INTRODUCTION 

The  main  parameters  of  seismic  signals  are  related  to  both  earthquake  sources  and 
propagation  medium  :  arrival  times  to  location  and  velocity  structure,  amplitude  to  source 
energy  and  attenuation  within  the  earth. 

If,  as  soon  as  the  beginning  of  the  century,  the  study  of  seismic  travel  time  is  tackled  in 
a  theoretical  way  by  Herglotz  and  Wiechert  (1909)  to  obtain  the  variation  of  the  velocity  with 
the  depth  from  readings  of  seismograms  and  by  Geiger  (1910)  to  localize  earthquake,  the 
study  of  the  amplitude  was  for  a  long  time  neglected. 

It  is  only  since  Richter’s  work  in  1935,  that  numbers  (i.e.:  magnitudes)  are  assigned  to 
quantify  the  size  of  earthquakes.  The  magnitude  is  based  on  amplitude  measurements  on 
instrumental  recordings.  The  original  purpose  was  to  facilitate  cataloging  of  earthquakes 
without  depending  exclusively  on  felt  intensities  to  compare  one  event  to  another.  Richter’s 
work  introduced  the  so-called  local  magnitude  scale  based  on  the  trace  amplitude  recorded  on 
the  Wood-Anderson  seismograph.  This  amplitude  is  a  function  of  the  ’size’  of  the  earthquake 
(magnitude)  and  of  the  attenuation  which  is  only  represented  by  an  empirical  formula  relating 
amplitude  and  distance.  Gutenberg  and  Richter,  during  the  1940’s,  extended  the  empirictd 
tables,  (i)  to  use  observations  made  at  distant  stations  and  on  seismographs  of  other  than  the 
standart  type,  (ii)  to  cover  earthquakes  of  significant  focal  depth  and  (iii)  to  enable 
independent  magnitude  estimates  to  be  made  from  body  and  surface  waves. 

It  is  necessary  to  wait  for  the  1960’s  to  have  a  new  impulse  in  the  amplitude  studies, 
with  the  model  ization  of  the  far- field  displacement  spectrum  produced  by  a  propagating  fault 
(Ben-Menahem  et  al.  1963,  Haskell  1966,  Brune  1970,  Madariaga  1976)  depending  on  the 
source  parameters  (seismic  moment,  fault  length,  stress  drop). 

The  end  of  the  1980’s  shows  a  real  burst  of  the  number  of  attenuation  studies  at 
regional  distances  in  various  regions  of  world  (for  references  see  Herraiz  et  al.  1986) 
following  the  increase  of  numerical  seismic  signal  recordings. 

The  attenuation  can  be  defined  by  the  wave  amplitude  decay  that  occurs  when  a  wave 
propagates  through  a  real  medium.  This  decay  has  two  very  different  origins:  (i)  the 
geometrical  spreading  which  is  independent  of  the  frequency  in  numerous  cases,  and  (ii)  the 
anelastic  attenuation  which  acts  generally  as  a  low-pass  filter. 

The  anelastic  attenuation  is  usually  parameterized  by  the  quality  factor  Q.  The  observed 
Q  is  thought  to  be  a  combination  of  two  different  loss  mechanisms:  a  frequency-independent 
intrinsic  quality  factor  Q,  and  loss  due  to  scattering  process. 

The  intrinsic  Q  is  given  by  1  IQ,  =  -AE/(2ttE),  where  AE/E  is  the  fraction  of  energy  lost 
in  a  wave  cycle  due  to  imperfection  in  the  elasticity  of  the  material. 

The  effect  of  the  scattering  is  not  exactly  energy  loss,  but  rather  a  redistribution  of 
energy  in  the  medium.  The  amount  of  scattered  energy  is  very  dependent  on  the  ratio  between 
obstacle  and  wavelength  sizes,  and  the  resulting  effect  is  strongly  frequency  dependent  The 
scattering  can  be  also  parameterized  by  a  quality  factor  given  bv  l/@s  =  f  Igv  (Dainty  1981) 
where  g  is  the  turbidity  (i.e.  the  inverse  of  mean  free  path),  v  is  the  wave  velocity,  and  /  is 
the  frequency. 

The  attenuation  is  the  sum  of  these  two  contributions  and  the  observed  Q  can  be 

written  as 
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The  knowledge  of  the  attenuation  can  have  three  different  purposes:  (i)  the  computation 
of  the  source  parameters,  (ii)  the  estimation  of  the  tectonic  activity  of  the  area  and,  (iii) 
perhaps,  the  earthquake  prediction. 

-The  seismic  signis  emitted  by  the  earthquake  are  strongly  affected  by  the  attenuation 
and  then  the  adequate  knowledge  of  the  attenuation  of  high-frequency  seismic  energy  within 
the  earth  is  very  important  in  order  to  study  the  seismic  source  (seismic  moment,  source 
radius  and  stress  drop)  and  to  predict  ground  motion  from  theoretical  sources  with  application 
to  seismic  hazard  analysis. 

-It  is  now  generally  accepted  that  observed  Q  increases  with  frequencies  between  1  and 
30Hz  (Mitchell  1981).  This  frequency  dependence,  which  usually  is  expressed  by  the  relation 
Q(f)  =  Qofa,  is  often  found  to  be  stronger  with  increasing  tectonic  activity  and  is  often 
explained  as  being  related  to  the  decrease  of  homogeneity  in  crust  (Aki  1981).  On  a  very 
large  Chinese  data  set  Jin  and  Aki  (1988)  found  a  remarkably  linear  relation  between  the 
logarithm  of  the  coda  ^  maximum  earthquake  within  the  area  sampled  by  coda 
waves.  Then  the  knowledge  of  the  attenuation  could  be  used  to  evaluate  the  maximum 
earthquake  in  a  region  and  would  be  very  useful  for  siting  critical  facilities  such  as  nuclear 
power  plants. 

-Attenuation  is  a  good  indicator  of  the  heterogeneous  state  and  of  the  dissipative 
property  of  the  lithosphere;  and  it  is  expected  that  these  conditions  in  and  around  the  seismic 
source  volume  change  before  an  earthquake.  Then  it  is  reasonable  to  search  for  temporal 
variation  of  the  attenuation  and  probe  the  possibility  of  using  this  variation  for  earthquake 
prediction  (Gusev  et  al.,  1984,  Jin  et  al.,  1986).  But  some  authors  have  found  no  precursory 
change  in  attenuation  (Scherbaum  1985,  Huang  1988)  and  this  use  of  the  attenuation  is  not 
now  operational. 

In  this  paper,  after  a  quick  review  on  the  two  main  methods  to  evaluate  the  quality 
factor,  we  present,  as  examples,  a  number  of  studies  of  attenuation  in  France. 

1-7.2  METHODS 

•  The  most  common  way  to  estimate  the  anelastic  attenuation  in  the  earth  consists  in 
studying  the  decay  of  the  seismic  energy  of  a  wave  with  distance.  The  quality  factor  obtained 
by  this  method  will  be  written  as  Qw  where  w  represents  the  name  of  the  studied  wave  (i.e., 
Pg,  Lg,  etc.)  This  method  needs  several  stations  in  order  to  have  a  good  sampling  of  the 
cpicentral  distances. 

•  It  is  also  possible  to  measure  the  attenuation  with  only  one  station.  The  tail  of  a 
wavetrain,  which  contents  the  energy  recorded  after  the  direct  wave  is  named  its  coda.  Two 
points  of  view  are  necessary  to  explain  the  amplitudes  of  these  two  wave  types  (direct  and 
coda): 

-  For  direct  waves,  a  deterministic  point  of  view  associates  without  ambiguity  one  effect 
to  one  well-known  heterogeneity. 

-  On  the  other  hand,  a  stochastic  approach  considers  coda-waves  as  diffracted  primary 
waves  by  a  uniform  distribution  of  heterogeneities.  In  this  case  the  medium  can  be  described 
by  a  set  of  statistical  parameters  (e.g.:  mean  free  path,  density  of  scatterers)  which  depend  on 
the  physical  theory  used  to  explain  the  diffraction.  In  all  the  cases  the  energy  loss  of  a 
primary  wave  can  be  summarized  by  the  use  of  a  quality  factor  Qs  representing  all  the 
scattering  effects,  which  is  surimposed  to  the  intrinsic  quality  factor  Q{. 

In  a  given  region,  the  relative  importance  of  the  absorption  and  the  scattering 
determines  the  shape  of  seismograms  of  local  earthquakes.  High  (2,  and  strong  scattering  is 
associated  with  very  long  time  duration  seismograms.  This  is  the  case  for  lunar  seismograms. 
On  the  other  hand,  a  low  Qi  (i.e.,  strong  absorption)  diminishes  the  time  duration  on  the 
recordings.  Hence  the  coda  of  seismograms  is  affected  by  attenuation  and  therefore  the  coda 
recordings  can  be  used  to  evaluate  attenuation.  The  quality  factor  obtained  by  coda 
measurements  will  be  identified  as  Qc. 

7.2.1  Qw  EVALUATION 

As  the  source  and  path  effects  are  separated,  the  spectral  amplitude  of  the  wave  w 
recorded  at  station  i  produced  by  earthquake  j  is  : 

Awij(f  -A)  =  Swij(f)  Pwij(f  ,A) 
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where  A  denotes  the  epicentral  distance,  /  the  frequency,  SwiAf  )  the  source  contribution  of 
the  jth  earthquake  and  Pwi.if  A)  the  propagation  effect.  "Each  factor  can  be  modeled 
(Campillo  et  al.,  1985)  as: 

-  The  source  effect  : 


SwlJif)  =  nj(J)Rl] 


•  is  the  source  amplitude  spectrum  of  the  jth  earthquake.  Following 

Brune  (1970),  this  spectrum  is  characterized  by  3  parameters:  0(0)  the  low- 
frequency  level,  which  is  proportional  to  the  seismic  moment  M  0;  fc  the  comer 
frequency  related  to  the  fault  length;  and  m  the  power  of  the  high-frequency 
asymptote.  The  spectrum  is  described  by  the  relation: 


Q(0) 

J  c 


1-7.1 


•  RtJ  is  the  radiation  pattern  of  the  jth  earthquake  to  the  ith  station,  function 
of  the  source  geometry  for  a  given  wave  type. 

-  The  propagation  effect  : 


Pwij  if  -A)  =  Ew  (A)  AAW  (f  ,A)  STWI  if ) 


•  Ew( A)  is  the  geometrical  spreading.  Since  epicentral  distance  range  is  small 
(<10°)  this  factor  can  be  written  in  the  approximative  form 

Ew(  A)  =  A~7  1-7.2 

•  AAW  if  ,A)  is  the  anelastic  attenuation  factor.  It  can  be  written  as: 

nf  a 

AAw(f,  A)  =  c  Q-(f)v  1-7.3 

with  Qw  if )  -  Qwy\fa  the  quality  factor  and  v  the  group  velocity  of  tire  wave. 

•  STwi(f )  is  the  station  response  and  represents  the  particular  characteristics  of 
the  station  which  are  not  included  in  the  instrumental  response.  In  many 
attenuation  studies  this  factor  is  neglectc  :  but  our  experience  convinces  us  of 
its  great  importance.  This  factor  is  mainly  the  frequency  response  of  the 
superficial  structures  beneath  the  station.  Generally  these  structures  are  not 
known  with  sufficient  accuracy  and  then  this  factor  must  be  deduced  from 
observations.  In  this  case  it  is  necessary  to  assume  for  each  frequency  f: 

,v 

nsr  „,(/)=  l  1-7.4 

I  =  1 

where  N  is  the  number  of  stations. 


After  linearization  by  taking  logarithms,  the  spectral  amplitude  A  of  the  wave  is: 
log(A,„y  if  ))  =  log(5,v</  (f ))  ~  Y  log(A)  -  ~~  +  log iSTwl  if  ))  1-7.5 

This  equation  can  be  solved  by  an  iterative  process  using  the  least  squares. 

7.2.2  Qc  EVALUATION 


7 


Many  models  have  been  used  by  the  investigators  in  coda-wave  analysis.  A  cornerstone 
of  these  models  is  the  relation: 


W.  0  =  W)  C(f,  t) 


1-7.6 


It  relates  the  power  spectrum  d>  of  the  coda  for  frequency  /  at  lapse-time  t,  with  a 
function  'P  which  includes  only  source  parameters  and  a  function  C  which  represents  the 
effects  of  a  large  geographical  area  and  is  independent  of  distance  and  details  of  path 
connecting  source  to  station.  This  relation  is  based  on  experimental  evidences  (Aki  and 
Chouet,  1975). 

It  also  indicates  the  two  steps  which  must  be  accomplished  in  a  coda-wave  study: 

-  finding  a  mathematical  expression  (model)  for  both  vf'(/)  and  C(f ,  t), 

-  relating  <£>(/ ,  t )  and  the  coda  amplitudes. 

•  Among  the  different  models  which  have  been  suggested  to  explain  the  formation  of 
coda-waves  (e.g.:  surface-waves  scattering,  diffusion  model,  multiple  scattering,  etc.),  the 
"Single  Backscattering  Model"  (Aki  and  Chouet,  1975,  Rautian  and  Khalturin,  1978)  seems  to 
be  the  most  common.  It  describes  the  coda-waves  as  a  superposition  of  backscattering  body- 
waves  from  discrete  scatterers. 

Under  these  assumptions,  the  function  describing  the  medium  effects  can  be  written: 

C(f,t)  =  r2ye-2nftlQ‘  1-7.7 

The  first  term  of  this  function  is  the  geometrical  spreading  factor  and  the  exponential 
term  represents  both  intrinsic  and  scattering  attenuation  effects. 

•  The  second  step  is  the  relation  between  the  power  spectrum  <J>(/,r)  and  the 
amplitudes  Ac(f ,  t)  of  the  coda  for  a  given  frequency  /  and  a  grien  lapse-time  t.  As  for 
other  problems  of  time-frequency  analysis,  two  approaches  can  be  followed  that  proceeds  in 
the  time  and  frequency  domain  respectively. 

•  In  the  time  domain,  the  method  consists  of  an  evaluation  of  the  mean  amplitude  over 
a  little  moving  window  of  signal  filtered  around  a  given  frequency. 

•  In  the  frequency  domain,  the  amplitude  is  obtained  by  an  averaged  value  around  the 
given  frequency  /  of  the  Fourier  spectrum  of  a  moving  time-window  centered  at  time  t. 

In  both  cases  the  amplitudes  are  related  to  the  power  spectrum  by  the  relation: 

Ac(f ,  t)  =  <Wf7tW  1-7.8 

where  A/  is  the  bandwidth  of  the  filter  in  the  time-domain  method,  and  the  bandwidth  of 
frequency  averaging  in  the  frequency  method. 

Combining  relations  [1-7.6],  [1-7.7]  and  [1-7.8]  we  get  the  classical  formula: 

Ac(f ,  t)  =  <WJ)Kf  r'<  e~nf,‘Q‘ 

Then,  the  Q  -factor  can  be  obtained  by  a  least-squares  fit  of  the  logarithm  of  this  formula. 

1-7.3  DATA 

It  results  from  the  preceding  section  that  both  Qc  and  Qw  evaluations  can  be  carried  out 
after  spectral  computations  of  part  of  signals.  Though  some  studies  have  been  carried  out 
using  analog  signals,  numerical  records  are  more  useful  because  they  allow  the  investigators 
to  compute  the  amplitude  spectra  using  the  Fast  Fourier  Transform  algorithm. 

Concerning  the  data  set,  Qw  studies  imply  the  use  of  a  network  in  order  to  have  a  good 
sampling  of  distance  for  each  event.  On  the  other  hand,  though  a  coda-f2c  value  can  be 
obtained  for  each  station-earthquake  couple,  the  use  of  a  network  allows  us  to  get  a  statistical 
value  for  the  studied  region  including  stations  and  earthquakes. 

In  France,  to  study  the  attenuation,  we  use  the  waves  recorded  by  the  short-period 
vertical  seismic  stations  of  the  seismic  network  run  by  the  Laboratoire  de  Detection  et 
Geophysique  (L.D.G.)  of  the  French  Commissariat  h  VEnergie  atomique  (C.E.A.)(fig.  1-7.1). 
Each  station  is  mainly  composed  of  a  1  Hz  high-gain  vertical-component  seismometer. 
Numerical  data  are  digitized  with  a  50  Hz  sampling  frequency  and  telemetered  to  the 
recording  center  where  they  are  recorded  on  paper  and  magnetic  tape  for  digital  processing. 
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The  spatial  distribution  and  the  magnitude  range  of  the  western  Europe  seismicity  (see 
figure  1-6.10  of  this  book)  allow  us  to  have  accurate  evaluation  for  both  Qw  and  Q;.  For 
Qw  measurements,  the  area  is  controlled  by  the  station  repartition  ana  then  its  evaluation  will 
be  made  only  for  Central  France.  A  better  spatial  repartition  on  France  will  be  obtained  by  the 
coda  measurements. 

1-7.4  RESULTS 

7.4.1  Qw  STUDY 

For  intermediate  distances  (150-1 500km''  along  continental  propagation,  seismograms 
are  mostly  dominated  by  crustal  phases  and  especially  by  Lg,  as  tar  as  amplitudes  and 
durations  are  concerned  (fig.  1-7.2).  Interpreted  as  a  complex  composition  of  multi-reflected  S 
waves  within  the  crust  (Bouchon  1982),  the  l.g  wavetrain  is  particularly  well-adapted  for 
attenuation  measurement. 

We  present  here  the  main  results  obtained  by  Campiiio  et  al.  (1985)  on  the  attenuation 
in  die  crust  beneath  central  France.  In  this  work  the  equation  [1-7.5]  is  solved  for  a  set  of  18 
earthquakes  (fig.I-7.2)  recorded  by  22  stations  of  the  permanent  L.D.G.  network.  In  order  to 
increase  the  stability  in  the  inversion  process,  the  geometrical  spreading  factor  (equation  [I- 
7.1])  must  be  fixed.  Using  numerical  simulations  of  Z^-wave  propagation  in  an  elastic 
homogeneous  layered  crust,  where  only  geometrical  attenuation  is  considered,  we  have  found 
y  =  0.83  (Campiiio  et  al.  1984). 

7.4.1. 1  STATION  RESPONSE 

The  transfer  functions  obtained  are  not  absolute  but  relative  to  the  mean  of  die  station 
set  (equation  [1-7.4]).  They  are  essentially  representative  of  any  attenuation  or  magnification 
for  high  frequencies  by  factors  up  to  3  and  as  shown  by  figure  1-7.3  they  are  quite  the  same 
for  Pg  and  L.  waves.  They  are  also  rather  well  correlated  with  seismic  noise  spectrum 
variation  (fig.  1-7.3).  This  result  is  important  because  it  means,  that  for  a  quick  evaluation  of  a 
site  effect,  it  is  possible  to  compare  only  noise  spectrum  between  the  s;te  and  a  close 
reference  point. 

7.4. 1.2  QUALITY  FACTOR 

The  quality  factor  is  computed  for  three  parts  of  the  Lg  wavetrain  defined  by  the  group 
velocity  windows: 

3.1  km  Is  <  v  <  3.6  km  Is  Lgl  (maximum  amplitude  ox  Lg) 

2.6  km  Is  <  v  <3.1  km  Is  Lg2  (early  part  of  the  coda)  1-7.9 

2.3  km  Is  <  v  <  2.6  km  Is  Lg  3 

Assuming  the  same  geometrical  spreading  (y  =  0.83)  and  the  same  frequency 
dependence,  the  Qwq  obtained  for  these  three  parts  are  very  close  with  only  a  slight  increase 
from  Lg  ]  to  Lg  3  : 

QLgl=  290*f052 

QLg2  =  3 10*/ 052 

QLg  3  =  330V0'52 

These  relations  mean  that  the  Lg  waves  and  their  early  coda  have  the  same  apparent 
attenuation  and  confirm  the  similarity  of  the  Q  factors  computed  from  S  waves  and  from 
their  coda. 


7.4. 1.3  SOURCE  SPECTRUM 
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The  knowledge  of  the  spectral  energy  decay  and  of  the  station  responses  can  be  used  to 
compute  the  source  excitation  of  each  earthquake.  The  source  spectrum  presents  a  dispersion 
(Fig.  1-7  4)  probably  due  to  a  rough  modeling  of  propagation  which  does  not  take  into  account 
lateral  variations  of  the  crustal  structures.  Nevertheless,  the  spectra  pointed  out  the  main 
features  of  the  source  amplitude  spectrum  (equation  [1-7.1]),  and  can  be  used  to  evaluate  the 
seismic  source  parameters  (seirruc  moment,  fault  length  and  stress  drop).  Usually  this 
estimation  is  made  only  with  sufficiently  close  stations  to  avoid  attenuation  corrections,  or  for 
large  earthquake  with  measurements  on  surface  waves  recorded  at  distant  stations.  The  use  of 
Lg  wave  allows  to  extend  this  calculation  to  the  regional  distance  ranges.  The  three  parts  of 
the  Lg  wavetrain  lead  to  similar  source  spectra  (figure  1-7.4),  and  in  the  case  of  saturated 
records  the  early  parts  of  the  coda  can  be  used  to  compute  the  source  parameters. 

7.4.2  Qc  CTUDY 

This  section  presents  the  results  obtained  on  the  ftequency  dependence  of  Qc  measured 
from  the  coda  of  regional  earthquakes.  Its  high  frequency  behaviour  is  investigated  up  to  100 
Hz.  The  frequency  dependence  exhibits  clearly  regional  variations.  Qc  will  be  compared  with 
the  previous  Q[_g  determination. 

7.4.2. 1  Qc  as  a  function  of  frequency 

To  study  the  high  frequency  behaviour  of  the  Qc -factor  in  a  wide  frequency  band,  we 
use  a  broadband  recording  system  at  station  Lormes  (LOR)  of  the  French  L.D.G.  network. 
The  output  signal  from  the  seismograph  is  separated  into  3  channels  after  filtering  in  different 
frequency  bands  and  with  amplification  and  digitalization  rate  adapted  to  each  frequency  band. 
This  system  presents  a  dynamic  range  large  enough  to  study  Qc  in  a  wide  frequency  band 
ranging  from  1  to  100  Hz. 

A  typical  record  of  this  system  is  shown  on  the  figure  1-7.5.  It  is  an  mL=  3.6  earthquake 
which  occurred  in  the  southern  part  of  the  Rhinegrabcn  at  27 0  km  from  station  LCR.  On  the 
3  filtered  traces  one  can  se  the  4  classical  phases  of  a  regional  earthquake.  The  crustal 
wavetrains  Pg  and  Lg  are  clearly  seen  on  the  low  frequency  channel.  But,  for  frequencies 
greater  the.*  *♦ 0  Hz,  r  ly  the  refracted  waves  Pn  and  Sn  are  well  developped,  so  that  for  these 
frequencies  the  coda  is  the  coda  of  Sn  waves. 

Despite  this  observation,  the  frequency  variations  of  the  quality  factor  follow  a  very 
regular  power-law  function  (Fig.  1-7.6  ): 

Qc  =  133/°  81  1-7.10 

which  loads  to  a  low  value  for  the  quality  factor  at  1  Hz  Qc0.  This  result,  rather  different 
than  the  network  averageu  value  obtained  for  Qig,  suggests  a  more  complete  study  of  the 
regional  variations  of  coAa-Qc. 

7. 4. 2. 2  Regional  variations  of  Qc 

In  order  to  study  the  regional  variations  of  the  Qc  factor,  we  have  processed 
automatically  a  set  of  earthquakes  having  a  local  magnitude  ranging  from  3.8  to  5.0.  The 
earthquake  set  is  separated  into  three  parts  according  to  their  locations.  For  each  set  Qc  is 
computed  (fig.  1-7.7  for  eastern  part).  We  found  the  following  results: 


Qc 

=  (210  ±  40)/(06°  ±0-,0) 

Eastern  part 

1-7.11 

Qc 

=  (290  ±  50)/ (0  54  *  0  08) 

Central  part 

1-7.12 

Qc 

=  (320  ±  50)/(052  ±O06) 

Western  part 

1-7.13 

7.4.3  QLg  AND  Qc  COMPARISON 

We  have  obtained  estimations  of  S-wave  attenuation  within  the  crust  beneath  France 
with  two  different  methods  and  for  different  regions. 
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•  In  the  common  studied  zone  (i.e.:  Central  France),  the  results  are  very  similar  (see 
equations  1-7.9  and  1-7.12).  Another  similar  value  is  obtained  for  the  Western  part  of 
France(equation  1-7.13).  This  confirms  the  hypothesis  of  Aki  et  al.  [1986],  which  indicates 
that  Qc  and  QLg  agree  very  well  as  soon  as  QLg  >  150. 

These  values  suggest  that  the  central  part  of  France  is  situated  between  active  and  stable 
zones  (Figure  1-7.8).  Following  Jin’s  relation  [1988],  the  maximum  magnitude  corresponding 
to  Qc o  =  300  is  roughly  6. 

•  For  the  Eastern  part  of  France,  both  Qc 0  and  the  frequency  dependence  suggest  that 
this  region  is  more  tectonically  active  (see  equations  1-7.10  and  1-7.11).  The  maximum 
magnitude  estimation  is  6.7  (with  Qc  =  210)  or  7.5  (with  Qc  —  133)  which  is  reasonable 
according  to  the  size  of  the  greatest  earthquakes  suffered  in  this  region  (Basel,  1356  and 
Swabian  Jura,  191 1). 

In  summary  this  study  shows  the  importance  of  the  regionalized  measurement  of  the 
Q  -factor  and  suggests  a  more  detailed  approach  using  a  discrete  model  of  attenuation  and  a 
tomographic  inversion. 

7.4.4  TOMOGRAPHY  OF  THE  QUALITY  FACTOR 

The  attenuation  measurements  made  from  L„  amplitudes  are  suitable  for  tomogTaphic 
studies  because  the  mean  value  of  Q  deduced  from  one  record  concerns  a  slice  of  crust 
around  the  source-receiver  path.  The  tomographic  approach  is  more  difficult  to  apply  when 
considering  Qc  measurements  because  the  volume  of  lithosphere  concerned  is  large. 

We  have  computed  the  distribution  of  Qs  in  the  crust  from  a  large  set  of  Lg  records  in 
different  frequency  ranges  (Campillo  ;  nd  Piantet,  1990).  We  use  raw  data  from  430 
earthquakes.  Our  results  confirm  that  the  eastern  part  of  France,  namely  the  Alps  region,  is 
associated  with  low  apparent  (2 -values  (figure  1-7.8). 

On  the  other  hand,  we  found  in  the  western  part  of  France,  a  zone  where  a  stronger 
attenuation  is  associated  with  a  high  reflectivity  of  the  lower  crust,  as  shown  by  a  wide  angle 
reflection  profile. 

The  tomography  reveals  clearly  the  heterogeneity  of  the  distribution  of  attenuation  and 
shows  that  some  anomalies  are.  difficult  to  predict  from  the  analysis  of  surface  geological 
structures. 


1-7.5  CONCLUSION 

The  quality  factor  of  both  coda-waves  and  Lg-waves  has  been  studied  using  a  set  of 
regional  events  recorded  on  a  seismic  homogeneous  network. 

•  Using  coda-Q  method,  the  frequency  dependence  of  the  Q  -factor  has  been  clearly 
shown  in  a  wide  frequency  band  (1  to  100  Hz). 

•  The  two  methods  lead  to  very  similar  results  in  the  common  studied  zone  (see 
equations  1-7.9  and  1-7.12)  for  the  studied  frequency  band  (1  to  15  Hz).  This  is  consistent 
with  other  studies  (e.g. :  Singh  and  Herrmann,  1983),  and  suggests  that  both  QLg  and  coda-Q 
measure  the  total  Q ,  resulting  from  intrinsic  and  scattering  attenuations. 

•  Then  these  attenuation  functions  have  been  used  to  compute  source  spectra  which  are 
the  basic  elements  for  the  evaluation  of  physical  seismic  parameters  such  as  seismic  moment, 
stress  drop  and  fault  length.  This  has  been  done  by  an  inversion  technique  using  a  model 
including  source  spectrum  pattern  and  station  responses.  We  have  shown  that  this  last  factor 
plays  a  prominent  part  and  cannot  be  neglected. 

•  Concerning  the  interpietation  of  Q  -values,  their  correlation  with  the  tectonic  activity 
has  been  established  in  France.  This  confirms,  at  a  regional  scale,  the  result  pointed  out  by  Jin 
et  al.  (1986),  concerning  the  ability  of  Q  -factor  to  be  a  useful  parameter  to  characterize  the 
tectonic  activity. 

Funner  studies  might  investigate  the  possibility  of  temporal  changes  in  Q  -values.  If  this 
fact  is  proved,  Q  -factor  will  be  also  a  worthy  parameter  for  earthquake  prediction,  which  is 
one  of  the  most  important  goals  of  present-day  seismology. 
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FIGURE  CAPTIONS 

Figure  1-7.1:  Map  of  stations  (triangles)  and  earthquakes  used  for  the  attenuation 
estimations.  Circles  refer  to  the  QL„  evaluation  and  crosses  to  the  Qc  study.  The  star 
shows  the  position  of  the  event  used  in  the  wide  frequency-band  study. 

Figure  1-7.2:  Examples  of  regional  earthquake  seismograms.  The  epicenter  is  located  in 
the  Pyrenees  (France).  The  seismograms  are  delayed  on  the  Pg  arrival-time.  For  clarity, 
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only  4  stations  are  displayed  among  the  16  which  have  recorded  this  event. 

Figure  1-7.3:  Examples  of  station  responses.  The  responses  are  calculated  from  L 
waves  (solid  line),  P  waves  (thin  line)  and  from  noise  (dashed  line).  This  figure  show* 
the  importance  of  this  factor  (magnification  up  to  3).  Note  that  the  different  evaluations 
(Lg ,  Pg  and  noise)  lead  to  similar  results. 

Figure  1-7.4:  Source  excitation  spectra  obtained  from  Lgi,  Lg2  and  Lg 3  waves.  Each 
solid  line  represents  the  spectrum  at  each  station  after  propagation  and  station  response 
corrections.  Heavy  lines  figure  the  best  fit  and  the  deduced  comer  frequencies  are 
shown. 

Figure  1-7.5:  Example  of  seismic  vertical  signal  for  a  mL  =  3.6  earthquake  recorded 
with  the  temporarily  implemented  broad-band  system  in  station  LOR.  The  distance  of 
this  event  is  270  km.  Each  trace  is  the  same  signal  filtered  with  the  bandwidth  displayed 
on  the  left  of  each  trace. 

Figure  1-7.6:  The  left  part  of  the  figure  shows  the  enveloppes  of  the  coda  for  different 
frequencies  referenced  on  the  right.  The  lines  are  the  results  of  the  least-squares  fit  of 
the  studied  part  of  the  coda  which  yields  to  Qc  -values  also  referenced  on  the  right.  The 
flat  level  connected  to  each  curve  by  a  dotted  line  shows  the  noise  level  in  the 
corresponding  frequency  band.  The  right  part  of  the  figure  shows  the  frequency 
variations  of  the  Qc  values  and  the  corresponding  least-squares  fit  by  a  power  law- 
function. 

Figure  1-7.7:  The  top  part  of  the  figure  shows  the  coefficients  of  the  power-law  function 
Qc(f)  =  Qofa  resulting  from  the  least-sqarred  fit  of  the  data  set  represented  on  the  map 
below.  Each  trace  connects  a  station  (triangle)  to  an  earthquake  (dot). 

Figure  1-7.8:  Tomographic  map  of  QLg  at  3  Hz  obtained  by  a  global  inversion  of 
attenuation,  station  responses  and  source  spectra  from  a  430  earthquakes  set  recorded  on 
the  L.D.G.  network.  The  strong  variations  of  Q  at  a  regional  scale  is  clearly  shown,  and 
indicate  that  they  must  be  taken  into  account  for  strong  motion  decay  evaluation. 
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ATTENUATION  OF  CRUSTAL  WAVES  ACROSS  THE  ALPINE  RANGE. 
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Although  structural  boundaries  are  known  to  affect  the 
propagation  of  short  period  regional  phases  ,  our  understanding 
of  the  phenomenon  remains  poor.  Several  reasons  may  be  invoked  to 
explain  this  gap.  The  first  one  is  the  extreme  variability  of  the 
situations  encountered  in  the  Earth  and  our  poor  knowledge  of  the 
geometries  of  the  deep  crustal  structures.  A  second  reason  is 
the  small  number  of  detailed  seismological  studies  of  regional 
phases  across  a  structural  boundary,  particularly  with  regard  to 
the  quantitative  analysis  of  amplitude  variations.  Finally,  our 
modeling  capabilities  have  lagged  behind  the  medium  complexity 
revealed  by  geophysical  investigations  and  structural 
interpretations  . 

On  the  other  hand  several  zones  of  anomaly  in  continental 
areas  have  been  reported  for  Lg  waves  such  as  in  the  Himalyan 
Belt  (Ruzaikin  et  al.,  1972)  or  in  the  North  Sea  (Kennett  and 
Mykkeitveit,  1984).  As  we  will  show  in  this  study,  a  zone  of 


strong  weakening  of  Lg  also  exists  in  the  western  Alps. 


This 
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region,  comprising  part  of  France,  Italy  and  Switzerland  has  been 
extensively  studied  by  geologists  and  geophysicists  (foi  a  review 
see  Vialon  and  Thouvenot  ).  Following  numerous  refraction  and 
wide  angle  seismic  experiments,  a  deep  reflection  profile  was 
recently  conducted  through  the  range  (Ecors  Project).  The  deep 
structure  along  a  East-West  cross  section  is  therefore  known  in 
its  gross  features. 

Seismological  evidences 

The  existence  of  an  anomaly  of  propagation  of  the  crustal 
phase  Lg  in  the  Western  Alps  is  directly  revealed  by  the  analysis 
of  the  seismological  bulletin  of  the  french  LDG  network.  For  a 
set  of  earthquakes  in  this  region,  we  have  drawn  cn  Figure  la  the 
ray  paths  for  which  Pn  phases  were  read  and  picked  while  m 
Figure  lb  are  plotted  the  ra,r  paths  for  which  Lg  arrivals  were 
picked.  These  two  maps  allr  us  to  compare  grossly  the 
characteristics  of  propagation  in  the  upper  mantle  and  in  the 
crust.  They  suggest  that  a  strong  anomaly  of  crccr ja f ion  of  the 
crustal  phase  Lg  exists  in  the  region.  Cn  the  maps  is  also 
represented,  by  a  thick  dotted  line,  the  axis  cf  the  strong 
Bouguer  gravity  anomaly  coresponding  to  the  Ivrea  Zone.  This 
simple  comparison  suggests  that  the  same  deep  structure  is 
responsible  for  the  extinction  of  Lg  and  for  the  positive  gravity 
anomaly . 
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In  order  to  cor.ii^m  the  presence  of  a  zone  of  extinction  of 
Lg  in  the  Alps,  and  obtain  quantitative  information  about  this 
phenomenon,  we  selected  a  series  of  earthquakes  ensuring  in  the 
region  and  we  computed  the  density  of  spectral  amplitude  for  Pn 
and  Lg .  The  Pn  amplitude  is  measured  in  the  group  velocity 
window  8. 2-7. 9  km/sec  while  the  Lg  amplitude  corresponds  to  the 
window  3. 6-2. 9  km/sec.  Examples  of  records  used  in  this  study 
are  presented  in  Figure  2.  Figure  2a  shows  a  record  obtained  at 
station  LOR  (Central  France:  47 . 16N  3.51E)  for  an  earthquake  in 
the  region  of  Gap  (44.  ON  6.22E),  in  the  western  part  of  the 
Alpine  range.  The  path  does  not  cross  the  Alps.  The  Lg 
wavetrain  is  clearly  visible  and  presents  a  serie  of  peaks  much 
larger  than  the  coda  of  earlier  arrivals.  A  spectral  analysis  of 
this  record  indicates  that,  in  the  spectral  domain  Lg  dominates 
Pn  for  frequencies  up  to  8  Hz.  On  the  contrary,  we  present  in 
Figure  2b  a  record  from,  an  earthquake  near  Genova  (44.52N  9.51E) 
obtained  at  station  BGF  in  Central  France  (46.33N  02.50E).  In 
this  case  the  travel  path  crosses  the  zone  of  positive  Bouguer 
anomaly.  The  Lg  wavetrain  is  no  longer  visible  on  the 
seismogram.  In  this  case,  the  spectral  amplitude  of  the  signal 
in  the  Lg  group  velocity  window  is  about  the  same  as  the  one  in 
the  Pn  window  for  frequencies  below  3  Hz,  while  at  higher 
frequencies  the  Pn  amplitude  is  larger. 

We  selected  a  set  of  earthquakes  in  Switzerland,  Northern 
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Italy  and  Southeastern  France  recorded  at  short  period  stations 
of  the  LDG  (Laboratoire  de  Detection  Geophysique,  France)  and  IGG 
(Istituto  Geofisico  di  Genova,  Italy)  networks.  The  location  of 
these  events  is  presented  in  Figure  3.  For  each  record  we 
computed  the  ratio  of  Lg  to  Pn  spectral  amplitudes.  To  present 
the  results  in  a  convenient  way,  we  defined  a  regular  grid  with  a 
mesh  of  20  km  that  covers  the  region  of  interest.  We  simply 
affected  to  each  cell  a  value  equal  to  the  mean  value  of  the 
ratio  of  amplitude  Lg/Pn  computed  for  all  the  paths  that  cross 
the  cell.  This  operation  was  performed  in  different  frequency 
ranges.  The  results  obtained  at  2  Hz  and  3  Hz  are  presented  in 
Figure  4.  These  images  indicate  the  presence  of  the  anomaly  of 
propagation  east  of  the  Alpine  range.  We  cannot  resolve  the 

extension  to  the  east  of  the  zone  where  Lg  are  sharply  attenuated 
because  all  the  italian  stations  of  the  IGG  network  are  located 
west  or  south  of  the  positive  Eouguer  anomaly,  that  seems  to 
correspond  to  the  eastern  limit  of  "normal"  propagation  of  Lg . 
It  is  important  to  note  that  this  zone  of  anomaly  dees  nor 
correspond  neither  to  the  region  of  the  highest  topographies  nor 
to  the  one  of  the  deepest  Moho,  that  are  clearly  located  west  of 
the  region  characterized  by  the  anomaly  of  crustal  propagation. 
One  may  also  notice  the  existence  of  a  zone  of  strong  weakening 
of  amplitude  ratio  in  the  gulf  of  Genova.  This  is  probably  an 
effect  of  the  well  known  vanishing  of  Lg  in  oceanic  areas  (Press 
and  Ewing,  1952).  However  the  propagation  is  efficient  in  the 
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coastal  region,  suggesting  that  the  structures  responsible  for 
the  extinction  of  Lg  do  not  extent  in  the  Southern  part  of  the 
chain . 

At  a  frequency  of  2  Hz  the  value  of  the  mean  ratio  between 
Lg  and  Pn  amplitudes  is  larger  than  4  in  the  Western  Alps  and  in 
central  France,  but  decreases  to  less  than  C.3  in  or  east  of  the 
zone  of  positive  Bouguer  anomaly.  This  gives  a  rough  indication 
of  the  order  of  magnitude  of  the  sharp  decay  of  amplitude 
occuring  along  paths  crossing  the  range  in  the  East  West 
direction:  the  amplitudes  are  divided  by  more  than  10. 

A  simple  structural  model. 

In  order  to  be  able  to  model  the  propagation  of  Lg  through 
the  Alps  we  need  to  make  some  simplifying  assumptions.  We  shall 
only  consider  paths  in  the  East-West  direction  that  cross  the 
central  part  of  the  range.  In  this  region  the  structures  are 
elongated  mostly  in  the  North-South  direction  which  allows  us  to 
reduce  the  study  of  the  effects  of  these  structures  to  a 
two-dimensional  problem.  Following  the  model  of  lithospheric 
overthrusting  established  by  Menard  and  Thouvenot  (1984)  ,  we 
built  up  a  simplified  cross  section  of  the  Western  Alps  (Figure 
5)  that  includes  the  most  important  features  likely  to  severely 
affect  crustal  wave  propagation:  the  geometry  of  the  Moho  and  the 
presence  of  deep  sedimentary  basins  on  both  sides  of  the  range. 

To  compute  synthetic  seismograms,  we  used  the  discrete  wavenumber 
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boundary  integral  method  proposed  for  irregularly  layered  medium 
by  Bouchon  et  al .  (1989).  We  limited  our  computations  to  the 

case  of  SH  waves.  We  assumed  the  source  to  be  located  15  km 
beneath  the  Po  basin,  as  indicated  by  a  star  on  Figure  5.  The 
values  of  shear  wave  velocity  and  density  are  the  following:  2.8 
km/sec  and  2.8  in  the  sediments,  3.5  km/sec  and  3.1  in  the  crust 
and  4.7  km/sec  and  3.3  in  the  upper  mantle.  We  assumed  the 
quality  factor  to  be  100  in  the  sediments  and  400  elsewhere. 
Because  of  limitations  in  computation  time,  we  considered  only 
frequencies  lower  than  1  Hz .  We  calculated  time  series  of  80  sec 
duration  at  equally  spaced  receivers  located  at  epicentral 
distances  ranging  from  100  to  400  km.  To  get  a  view  of  what 
should  be  "normal"  propagation  we  also  present  in  Figure  6  a 
seismic  section  corresponding  to  the  case  of  a  flat  layered 
medium.  Both  sets  of  seismograms  are  plotted  with  a  reduction 
velocity  of  3.5  km/sec.  As  indicated  in  Figure  5,  point  (A) 
denotes  the  Western  limit  of  the  Po  sedimentary  nasin  while  point 
(B)  gives  the  location  of  the  eastern  edge  of  the  Rhone  and 
molassic  basin.  The  seismograms  computed  for  the  Alps  model  show 
a  spectacular  decrease  of  Lg  amplitude  in  the  central  part  of  the 
mountain  range.  However,  the  amplitude  strongly  increases  when 
the  waves  reach  the  molassic  sedimentary  basin  to  the  west  of  the 
Alps.  Two  causes  can  account  for  this  result.  First  the 
curvature  of  the  Moho  in  its  deepest  part  will  tend  to  focus  the 
reflecting  shear  waves  which  make  up  the  Lg  wavetrain.  Secondly, 
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the  diffraction  of  the  crustal  waves  near  the  edge  of  the 
molassic  basin  will  result  in  the  excitation  of  waves  trapped  in 
the  sediments  that  are  associated  with  high  amplitudes  as  shown 
in  Campillo  (1987).  The  importance  of  this  last  phenomenon  in 
our  simulation  is  shown  by  the  apparent  velocitv  of  the  late 
arrivals.  In  the  case  of  a  flat  layered  model,  the  arrivals  at 
locations  farther  away  than  (B)  have  apparent  velocities  higher 
than  the  reduction  velocity  of  3.5  km/sec,  in  accordance  with 
their  nature  of  multiply  reflected  waves.  For  the  Alpine  crustal 
model,  the  arrivals  in  the  same  distance  range  are  characterized 
by  apparent  velocities  lower  than  the  reduction  velocity, 
indicating  that  they  correspond  to  waves  trapped  in  the  low 
velocity  sediments.  The  increase  in  amplitude  and  duration  of 
the  wavetrain  when  penetrating  in  the  zone  of  sediments  may  be 
due  to  the  well  known  effect  of  amplification  by  a  soft 
superficial  layer.  Because  the  attenuation  in  sediments  is 
large,  one  may  argue  that  these  waves  will  rapidly  decrease. 

In  order  to  identify  the  part  corresponding  to  waves  that  sample 
the  entire  crust  (the  Lg  waves),  we  repeated  the  computation 
using  similar  models  where  the  sedimentary  layers  were  replaced 
with  crustal  materials.  The  seismograms  obtained  are  presented 
in  Figure  7.  They  indicate  that  multiply  reflected  waves 
continue  to  exist  after  the  crossing  of  the  range,  in  spite  of 
the  very  complex  geometry  of  the  crust-mantle  boundary.  This 
result  means  that  the  structural  model  proposed  does  not  produce 
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a  complete  extinction  of  Lg. 

To  measure  the  effective  decay  of  amplitude  due  to  the 
crossing  of  our  simplified  Alps  model,  we  present  in  Figure  8, 
the  spectral  amplitude  at  frequencies  around  0.75  Hz  as  a 
function  of  distance.  The  curves  are  plotted  both  for  the  flat 
model  and  for  the  model  of  the  Alps  including  the  sedimentary 
basins.  The  comparison  between  the  two  curves  indicates  that  the 
complex  geometry  of  the  structure  is  responsible  for  a  decrease 
by  a  factor  of  at  least  2  of  the  amplitude,  in  spite  of  the 
presence  of  the  sediments.  On  the  other  hand,  this  effect  is  not 
sufficient  to  account  for  the  observed  extinction  of  Lg. 

If  one  accepts  the  fact  that  the  main  features  of  the 
crustal  structures  beneath  the  Western  alps  are  known,  the 
existence  of  a  zone  of  extinction  of  Lg  can  be  explained  only 
partially  by  the  effect  of  large  scale  lateral  structural 
variations.  To  explain  the  observed  extinction,  we  need  to 
invoke  causes  which  were  not  taken  into  account  in  our  simplified 
model.  Obviously,  the  existence  in  the  crust  of  a  zone  of 
intense  attenuation  of  body  waves,  due  either  to  anelasticity  or 
to  scattering  on  small  scale  inhomogeneities,  could  be 
sufficient.  Quantitatively,  the  attenuation  required  corresponds 
to  a  100  km  wide  zone  with  a  mean  Q  of  100. 

Another  simplication  introduced  in  our  model  is  the 
assumption  that  the  crust  mantle  boundary  is  smooth.  The 
potential  effect  of  the  roughness  cannot  be  addressed  at  this 
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stage  since  we  are  not  able  to  quantify  the  most  likely  cause  of 
vanishing  i.e.  the  intrinsic  attenuation  of  the  crustal 
materials . 
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Figure  caption 


Figure  1 

a)  Ray  paths  drawn  for  Pn  waves  from  the  epicenter 
location  (*)  to  the  seismic  station  (  /\  1 

b)  Ray  paths  for  Lg  waves 

These  two  maps  allow  us  to  compare  grossly  the  characteristics 
of  propagation  in  the  upper  mantle  and  in  the  crust. 

Figure  2 


a)  Record  obtained  at  LOR  (Central  France:  H  7°,  16  N 
3°  51  E)  for  an  earthquake  in  the  region  of  Gap  in  the  western 
part  of  the  Alpine  range.  The  path  does  not  cross  the  Alps.  The 
Lg  wave  train  is  clearly  visible  and  presents  a  serie  of  peaks 
much  larger  than  the  coda  of  earlier  arrivals. 

b)  Record  from  an  earthquake  near  Genova  obtained  at 
BGF  in  Central  France  (H  6*  33N  02*50E).  In  the  case  the  travel 
path  crosses  the  zone  of  positive  Bouguer  anomaly.  The  Lg 
wavetrain  is  no  longer  visible. 

Figure  3 

Earthquakes  and  seismic  stations  (  /S  ) 

Figure  4 

Ratio  Lg/Pn  amplitudes  in  each  block: 

a)  at  f=  2  Hz 

b)  at  f=  3  Hz 

These  images  indicate  the  presence  of  an  anomaly  of  propagation 
east  of  the  Alpine  range.  Notice  the  existence  of  a  zone  of 
strong  weakening  of  amplitude  ratio  in  the  gulf  of  Genova. 

Figure  5 


Simplified  cross  section  of  the  western  Alps  used  in 
the  seismograms  modelization  with  deep  sedimentary  basins  on 
both  sides  of  the  range.  We  assume  the  source  (*)  located  15  km 
beneath  the  Po  basin. 

Figure  6 


a)  Seismograms  computed  with  a  "normal  propagation" 
in  a  flat  layered  medium. 
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b)  Seismograms  computed  with  the  "Alps  model"  showing 
a  spectacular  decrease  of  Lg  amplitude  in  the  central  part  of 
the  mountain  range. 

Figure  7 


Computation  using  similar  model  as  in  Fig. 6  where 
sedimentary  layers  are  replaced  with  crustal  materials.  The 
sismograms  indicate  that  multiply  reflected  waves  continue  to 
exist  after  the  crossing  of  the  range. 

Figure  8 

Effective  decay  of  amplitude  due  to  the  crossing  of 
the  Alps  models=  spectral  amplitude  at  0,75  Hz  as  a  function  of 
distance. 
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Objectives 

The  objectives  of  this  work  are  to  further  our  understanding  of  the  source  region 
geological  structure  effects  on  the  teleseismic  displacement  recorded  by  a  thin  aperture 
network.  The  calculations  are  carried  out  using  the  discrete  wavenumber  -  boundary 
integral  equation  formulation  linked  to  the  reciprocity  theorem. 

Research  accomplished 

We  investigate  here  the  source  region  effects  that  must  be  taken  into  account  when 
analysing  the  Yucca  Flat  nuclear  detonations  recorded  at  distances  of  80°  by  a  seismic 
network.  We  try  to  separate  the  relative  effects  of  the  structure  heterogeneities  and  of 
the  source  medium  which  contribute  to  these  teleseismic  ground  motions. 

The  farfield  displacement  is  assumed  to  result  from  the  convolution  of  the  source 
function,  i.e.  s(x),  by  the  propagation  path  response,  i.e.  p(X),  by  the  receiver  and  source 
site  functions,  i.e.  c(x),  and  by  the  recording  system  response,  i.e.  a(x).  Thus 

r(t)  =  s(x)  *  p(x)  *  c(x )  *  a(x) 

Different  formulations  of  the  source  function  s(t)  have  been  derivated  in  term  of  (i)  the 
reduced  displacement  potential  by  Haskell  (1967),  von  Seggern  and  Blandford  (1972), 
and  Helmberger  and  Hadley  (1981).  (ii)  Another  source  formulation  has  been  given  by 
Mueller  and  Murphy  (1971).  This  model  takes  into  account  both  the  yield  and  the 
detonation  depth.  The  model  scaling  factors  depending  on  the  medium  of  burial  have 
been  empirically  determined  by  Murphy  (1977)  foi  different  source  environments  or 
have  been  computed  by  Power  (1974)  using  a  code  that  simulates  underground  nuclear 
explosions  theory  for  tuff  having  various  porosities  and  liquid  densities. 

Different  numerical  methods  have  been  investigated  to  study  the  influences  of  the  near 
source  site  heterogeneities  on  the  teleseismic  P-wave  seismograms.  If  the  lateral 
heterogeneities  in  the  source  region  have  effects  at  teleseismic  distances  on  the  mb 
estimates  that  can  be  averaged  or  canceled  indeed  by  a  network  of  stations  at  multiple 
azimuths,  these  source  site  structures  are  of  great  importance  for  thin  aperture 
observations  (McLaughlin  et  al.,  1987,  McLaughlin  and  Jill,  1988).  Hasegawa  (1971  and 
1972)  using  a  composite  method  connects  the  recorded  seismograms  shape  to  the 
geological  complexity  of  the  source  environment.  Bouchon  (1976)  links  the  Thomson- 
Haskell  method  (Thomson,  1950;  Haskell,  1953)  with  the  reciprocity  theorem  and  shows 
the  great  influence  of  the  source  depth  and  of  the  medium  of  burial  on  the  mb  estimates 
for  events  in  layered  media.  Aki-Larner  (1970)  calculations  for  Yucca  Flat  explosions 
(Ferguson,  1988),  and  other  data  analysis  methods  (Taylor,  1983;  Lay  and  Wlec,  1987; 
Lay,  1987a  and  1987b)  have  been  processed  and  emphasize  the  near  source  structure 
effects  on  the  farfield  displacement  records. 

The  near  source  structure  effects  are  hereafter  computed  using  the  discrete  wavenumber 
-  boundary  integral  equation  formulation  associated  to  the  reciprocity  theorem.  This 
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discrete  wavenumber  method  has  been  previously  described  and  developped  in  the  P, 
SV,  and  SH  cases  to  study  the  effect  of  laterally  varying  media  onto  the  field  diffracted 
at  a  few  tens  of  kilometers  out  of  the  source  (Bouchon  and  Aki,  1977;  Bouchon  et  al., 
1989;  Gaffet  and  Bouchon  1989  and  1990).  We  study  the  farfield  source  environment 
influence  and  therefore  ignore  the  receiver  structures  for  such  a  problem  has  been  the 
subject  of  a  great  deal  of  works  which  give  prominence  and  describe  using  observed 
teleseismic  recordings  the  wave  diffraction  effects  that  are  due  to  the  topography  (Key, 
1967;  Bouchon,  1973;  Sanchez-Sesma  and  Rosenblueth,  1979;  Bard,  1982;  Kawase  and 
Aki  1988)  or  to  the  underground  structures  (Bard  and  Bouchon,  1980  and  1985)  located 
near  the  seismometer. 

A  profile  oriented  in  the  azimuthal  direction  of  France  has  been  infered  from  the 
Paleozoic  basement  depths  obtained  by  Taylor  (1983),  Ferguson  (1988)  and  McLaughlin 
et  al.  (1987)  for  the  Yucca  Flat  valley.  In  the  aim  to  determine  the  influence  of  the 
different  layers,  the  basin  was  firstly  made  up  of  a  wet  tuff  rocks  filling  (say  model  a), 
and  then  with  a  wet  tuff  volcanic  rocks  underlying  a  dry  alluvium  filling  (say  model  b). 
A  third  model  (say  model  c)  should  take  into  account  a  third  level  ( i.e .  the  dry  tuff 
level)  but  has  not  been  yet  processed  for  long  duration  computation.  The  model  b  and 
the  elastic  parameters  are  depicted  figure  1.  The  crosses  indicate  the  true  depths  of  the 
two  interfaces  considered.  The  horizontal  scale  gives  the  distance  in  kilometer  from  an 
arbitrary  origin.  II  detonations  are  studied  and  described  in  the  table  below.  The  map 
figure  1  displays  their  repartition  nearly  SW-NE  along  the  profile  given  figure  1. 


No 

Event  names 

Date 

Time 

'  N 

°  W 

Depth 

m* 

1 

PYRAMID 

16  Apr  1980 

20  00  00.1 

37  06  04 

116  01  50 

... 

579  m 

53 

2 

JORNADA 

28  Jaii  1982 

16  00  00.1 

37  05  29 

116  03  04 

640  m 

1 

3 

turquoise: 

14  Apr  1983 

19  05  00.1 

37  04  22 

116  02  45 

533  m 

KIM 

n 

MUNDO 

01  May  1984 

19  05  00.1 

37  06  22 

116  01  21 

567  m 

53 

5 

CAPROCK 

31  May  1984 

13  04  00.1 

37  06  1 1 

116  02  53 

600  m 

5.8 

6 

BRETON 

13  Sep  1984 

14  00  00.0 

37  05  12 

116  04  16 

483  m 

5.0 

a 

HERMOSA 

02  Apr  1985 

20  00  00. 1 

37  05  4 1 

116  01  56 

640  m 

n 

8 

GLENCOE 

22  Mar  1986 

16  15  00.1 

37  04  59 

116  03  58 

600  m 

5.1 

9 

TAJO 

05  Jan  1986 

15  04  00.1 

37  05  54 

116  00  56 

500  m 

53 

10 

GASCON 

14  Nov  1986 

16  00  00  1 

37  06  02 

116  02  53 

600  m 

5.8 

II 

BULLFROG 

30  Aug  1988 

18  00  00.1 

37  05  09 

116  04  07 

500  m 

5.0 

The  explosions  have  been  projected  onto  the  profile.  The  layout  shows  4  events  at  the 
same  horizontal  location  and  different  depths  (i.e.  Breton  (483  m),  Bullfrog  ( 500  m), 
Turquoise  (533  m),  and  Glencoe  (600  m))  and  a  set  of  7  events  from  west  to  east  (resp. 
Jornada,  Gascon,  Caprock,  Hennosa,  Pyramid,  Tajo,  and  Mundo).  The  mean  observing 
distances  is  80°  from  Yucca  Flat  and  corresponds  to  an  angle  of  incidence  from  the 
source  region  of  16.95°  (Pho  and  Behe,  1972).  Synthetic  seismograms  have  been 
computed  over  a  10  s  duration  window  at  a  such  distance  for  the  simple  basin  a  and  for 
the  two  layered  basin  b.  The  results  are  respectively  presented  on  the  left  and  right 
columns  of  the  figure  2.  The  farfield  basin  impulsive  response  has  been  convolved  with 
the  non  frequency  dispersive  Ricker  function  at  /  Hz.  The  time  domain  maximum 
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displacement  amplitude  increases  with  the  depth  of  burial  for  both  models  (See  Breton, 
Bullfrog,  Turquoise,  and  Glencoe).  The  general  feature  observed  is  the  amplitude 
increase  from  west  to  east  up  to  Pyramid  in  connection  with  the  depth  of  the  Palezoic 
basement.  This  shape  effect  is  magnified  by  the  presence  of  the  alluvium  layer  for 
which  the  maximum  displacement  location  always  remains  at  Pyramid  location.  The 
seismograms  are  all  the  more  broadened  because  the  alluvium  deposit  thickens.  Tajo  and 
Mundo  do  not  show  significant  variations  between  the  two  models  for  the  too  thin 
alluvium  coverage  above  them. 

We  have  computed  the  seismograms  that  should  be  recorded  in  France  for  the  11 
assumed  fully  coupled  detonations  and  for  yields  of  20  kt  and  120  kt.  The  Mueller- 
Murphy  source  (1971)  is  used  with  the  parameters  corresponding  to  the  wet  tuff-rhyolite 
medium  of  burial  (Murphy,  1977).  The  /*  value  is  0.7  (Cormier,  1982).  The  frequency 
responses  are  presented  figure  3  and  the  resulting  signals  ( i.e .  s( x)  *  p(x)  *  a( x))  are 
displayed  in  the  two  bottom  boxes  for  an  infinite  space  having  the  characteristics  of  the 
wet  tuff  layer.  The  main  frequencies  are  1  Hz  and  0.7  Hz  for  respectively  the  20  kt  and 
the  120  kt  shots  at  500  m  depth.  Besides  the  computed  seismograms  of  Tajo  and 
Mundo  which  are  scarcely  modified  in  comparison  with  the  half  space  solution,  the 
common  feature  is  an  enlargement  of  the  pP- wave  (figure  4).  For  both  yields,  the 
amplitude  increases  while  the  depth  of  the  Paleozoic  basement  decreases,  the  maximum 
amplitude  being  recorded  in  the  camber  region  in  the  eastern  part  of  the  valley.  If  the 
amplitude  seems  to  be  correlative  to  the  external  shape  of  Yucca  Flat,  it  appears  that  the 
frequency  of  the  maximum  amplitude  has  to  be  related  to  the  thickness  of  the  alluvium. 
The  figure  5  sums  up  these  results.  This  display  depicts  the  Amb  predicted  by  the 

synthetic  data.  The  A rnb  values  correspond  to  Log](j(A„  k  I  T)  -  |  Log]0  (/t,„ „  /  T)  j  , 

i.e.  the  difference  between  the  magnitude  of  each  detonation  k  with  a  reference  event  R 
which  predicts  the  maximum  magnitude.  The  general  behaviour  of  the  maximum 
amplitude  is  similar  for  the  two  yields  with  a  Arnb  variation  well  correlated  to  the 
Paleozoic  basement  shape  and  a  regular  decrease  of  the  period  of  the  maximum 
amplitude  with  the  thinning  of  the  alluvium.  A  source  effect  is  clearly  shown  for  the 
explosions  6  (Breton),  11  (Bullfrog),  3  (Turquoise),  and  8  ( Glencoe )  for  which  the  Amb 
behaviour  is  reversed  when  the  yield  increases.  Thus  for  small  yields  (e.g.  around  20  kt) 
the  displacement  amplitude  decreases  with  the  depth  of  burial  while  it  increases  for  big 
yields  (e.g.  around  120  kt). 

We  can  apply  these  synthetic  A mb  variations  to  estimate  the  correlative  yield  miss- 
estimation.  Different  formulas  may  be  used  to  estimate  the  yield  in  the  Yucca  Flat  test 
region,  i.e. 

mb  =  3.79  +  0.85  Log  (0  W  (Murphy,  1977)  (1) 

mb  =  3.71  +  0.89  Log  10  W  (Bache,  1982)  (2) 

These  two  relations  can  be  summarized  by  mb  =  a  +  b  Log^W.  Thus  on  the  one 
hand,  if  the  observed  magnitude  rnb  obs  is  smaller  than  the  magnitude  mb  lrue 
corresponding  to  the  true  value  of  the  yield  Wlrue,  the  miss-estimated  percentage  of  the 
yield  is 

- 1  A/n»  I 

I  AW  I  /  W,rue  =  1-10  b 

For  instance,  with  A mb  -  -0.25  the  underestimated  yield  is  half  times  the  true  yield.  On 
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the  other  hand,  if  the  observed  magnitude  mb  obs  is  greater  than  the  magnitude  mb  lrue, 
the  miss-estimated  percentage  of  the  yield  is  given  by 

I A mb  I 

l  AW  I  /  Wtrue  =  10  b  -  1 

Thus  with  A mb  =  +0.25,  the  estimated  yield  is  almost  2  times  the  true  yield.  The  b 
values  proposed  in  the  relations  1  and  2  are  close  together  so  that  the  predictive  yield 
estimation  effects  are  similar  for  both  magnitude  va\  energy  relations. 

Conclusion 

We  have  implemented  the  discrete  wavenumber  -  boundary  integral  equation  method  to 
compute  the  P- wave  farfield  response  of  structure  heterogeneities  in  the  source  region. 
The  computation  have  been  processed  using  the  reciprocity  theorem  to  study  11 
underground  explosions  occured  in  the  central  region  of  the  Yucca  Flat  valley,  Nevada, 
and  recorded  in  France.  We  have  shown  the  strong  influence  of  the  nature  and  of  the 
shape  of  the  basin  filling.  On  the  one  hand,  the  maximum  displacement  amplitudes  are 
connected  together  with  the  shape  of  the  valley  while  on  the  other  hand  the 
corresponding  frequency  is  correlative  to  the  thickness  of  the  alluvium  deposit.  The  mb 
anomalies  determined  in  this  work  would  be  used  to  reduced  the  uncertainty  in  the  yield 
estimations  made  with  data  obtained  from  a  thin  aperture  network.  The  complete  model 
which  has  not  yet  been  implemented  ( i.e .  the  model  c)  should  increase  the  model  b 
amplitude  behaviour  and  should  lengthen  the  first  arrivals  recorded  for  the  events  that 
occur  in  the  eastern  part  of  the  valley  where  the  dry  tuff  level  is  the  thickest. 
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11  in  the  tex.  are  depicted  by  the  black  circles  in  the  tuff  level. 
On  the  right  side  is  displayed  a  simplified  map  of  the  studied 
events.  The  gray  ellipse  is  plotted  to  emphasize  the  global  distri¬ 
bution  of  the  events  roughly  turned  toward  the  France. 
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Figure  2  -  Displays  of  the  synthetic  modellings  obtained  for  the  model  a  (left)  and  b 
(right)  for  the  /  /  studied  detonations.  The  source  is  for  both  eases  a  Ricker  having  a  l 
Hz  central  frequency  in  displacement.  The  relative  amplitudes  arc  written  in  the  left  side 

of  each  seismogram. 
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Figure  3  -  The  shapes  of  the  Mucller-Murphy  (1971)  source  that  should  be  recorded  at 
lelcscismic  dismnee  in  an  infinite  homogeneous  space  are  given  for  the  yields  20  kt  and 
120  kt  in  the  two  bottom  boxes.  In  these  boxes,  die  middle  vertical  thin  line  shows  the 
arrival  of  the  7  Hz  frequency.  The  transfcri  functions  of  the  Mucller-Murphy  (1971) 
source  function  (i.e.  modulus  IMM(/)I,  phase  d^MMO1))),  of  the  Futtcrman  (1962)  manlic 
attenuation  (i.e.  modulus  IFu(/)l),  of  the  recording  system  (i.e.  modulus  \S(f)\,  phase 
<t>[S(/))),  and  of  the  far  field  displacement  recording  (i.e.  modulus  ISis(/)i)  arc  depicted  in 
the  top  part  of  the  display. 
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Figure  4  -  Seismograms  obiainctl  for  the  Mucllcr-Murphy  source  (1971)  and  for  the  two 
yields  20  kt  (left)  and  120  kt  (right).  The  relative  amplitudes  and  their  corresponding  fre¬ 
quencies  are  written  at  the  left  of  each  seismogram.  The  seismograms  are  the  ones  calculat¬ 
ed  for  the  model  b.  The  solutions  obtained  in  a  half  space  arc  given  in  slight  superimposed 
dotted  lines. 
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Figure  5  -  Representation  of  the  synthetic  A mb  predictions  deduced  from  the  cal¬ 
culations  presented  figure  4.  The  geometrical  configuration  is  given  in  the  top  of 
the  figure.  The  6mb  (blade  squares)  and  the  frequencies  (white  circles)  at  which 
these  variations  have  been  computed  are  plotted  for  die  yields  20  kt  and  120  kt. 
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IV 

-  STUDY  OF  SEISMIC  SWARM  FROM  EASTERN  KAZAKH 

INTRODUCTION 

In  theory  a  seismic  event  can  be  defined  completely  by  a  time  function  including  : 

-  focal  mechanism 

-  magnitude 

-  location 

The  only  available  information  to  define  these  parameters  is  the  recorded  signal  at  a  station. 
This  signal  is  modified  along  its  path,  so  in  addition  to  the  information  about  the  source 
characteristics,  the  signal  provides  information  about  the  earth  structure  through  which  it 
travelled. 

In  this  study,  the  data  are  teleseisms  from  Eastern  Kazakh  recorded  by  the  french  network. 


STUDIED  DATA 


The  records  studied  here  were  provided  by  the  L.D.G.  (Laboratoire  de  Geophysique).  The 
french  seismic  network  consists  of  27  stations  with  the  same  features  :  the  velocity  response  is 
flat  above  1.0  Hz,  and  the  fundamental  period  of  the  instrument  is  1.0  Hz. 

Nineteen  events  from  Eastern  Kazakh  between  1981  and  1984  were  studied.  For  all  of  them, 
we  obtained  data  from  both  L.D.G.  and  I.S.C.  about  the  magnitude  and  the  location  (see 
table  1). 

The  locations  from  the  two  networks  were  somewhat  different  (see  figure  1)  ;  a  systematic  bias 
appeared  in  the  L.D.G.  locations,  probably  due  to  the  limited  range  of  azimuth  of  the  events. 
The  magnitudes  also  were  compared  and,  with  respect  to  the  I.S.C.  magnitude.  The  L.D.G. 
magnitudes  for  event  magnitudes  larger  than  6.0  are  over-estimated  whereas  they  are  under¬ 
estimated  for  smaller  magnitudes. 

In  this  study  I.S.C.  locations  and  magnitudes  were  used  for  reference. 
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MAGNITUDE  ESTIMATION 

One  of  our  goals  was  to  provide  a  way  to  estimate  magnitude  based  on  a  more 
"reliable"  method,  using  only  the  digital  signal  records  instead  of  using  graphic 
records.  We  first  supposed  that  a  recorded  signal  can  be  represented  as  the 
sum  of  K  times  an  unsealed  source  function  plus  noise  uncot related  with  the  signal 
and  characteristic  of  the  station. 

To  apply  this  relation,  some  assumptions  are  necessary: 

-  all  the  courses  have  to  be  the  same  (same  focal  mechanism)  and  in  our  case  we 
have  only  underground  nuclear  explosions  coming  from  the  same  region; 

-  the  K  factor  is  characteristic  for  each  station; 

-  the  aperture  of  the  network  has  to  be  small  with  respect  to  the  epicentral 
distance. 

Then,  the  K  factor  is  equivalent  to  A/T  in  the  theoretical  formula  for  the  magnitude 
and  could  be  defined  as  the  largest  eigenvalue  of  the  cross-correlation  matrix. 

Because  of  the  level  of  noise,  we  filtered  the  data  between  1.0  and  8.0  Hz.  T<~ 
position  the  signals  on  the  first  arrival,  we  used  a  cross-correlation  program. 

A  very  important  feature  of  the  signals  is  that  they  are  very  similar  for  different 
events  recorded  at  the  same  station  (see  Figure  2)  anu  for  this  case  the  cross¬ 
correlation  program  gives  good  results.  But  they  are  very  different  for  the  same 
event  recorded  at  different  stations.  This  result  shows  the  strong  influence  of 
the  receiver  function. 

The  construction  of  the  cross-correlation  matrix  requires  that  the  set  of  signals  used 
in  the  construction  be  complete.  That  is,  the  set  must  include  one  signal  for  each 
event-station  pair  in  the  set. 

Two  sets  of  data  were  used:  one  with  all  the  stations  (8)  having  recorded  all 

nineteen  events,  and  a  second  set  havi  tg  only  the  events  (8)  recorded  at  all  27  stations. 

With  these  calculations,  we  were  able  to  compute  a  K  factor  which  was  representative 
of  the  yield  for  a  given  station. 

Because  we  assumed  that  this  factor  was  equivalent  to  A/T,  we  cere  able  to  estimate 
a  magnitude  for  all  the  events  from  the  relation: 

mb  =  a  log(K)  +  b 

where  'a'  and  'b'  were  computed  by  least  squares  using  ISC  magnitude  for  reference. 

In  the  theoretical  formula  for  the  magnitude: 

mb  =  a  log(A/T)  +  b 

'a'  is  supposed  to  be  equal  to  1.0.  In  our  calculation  of  both  values  'a'  and  'b' 
we  found  that  'a'  value  varied  in  a  large  range  (from  0.77  to  1.23).  Since  we 
computed  the  values  for  all  the  stations,  we  were  able  to  map  'a'  values  as  a  function 
of  the  station  location. 

As  shown  or.  Figure  3,  the  correlation  between  the  location  and  'a'  value  is  strong, 
for  example  the  stations  located  in  the  South  of  France  show  the  highest  values  for 
'a'  (above  1.0)  while  the  stations  in  Normandie  show  a  value  around  0.88.  For  the 
stations  in  Massif  Central  and  Morvan,  the  relation  is  less  obvious  because  of  the 
geologic  complexity. 

Thus,  it  seems  important  for  accurate  computation  of  magnitude  to  take  into  account 
the  fact  that  the  'a'  value  could  be  very  different  from  1.0  and  is  strongly  influenced 
by  the  station  environment. 

With  all  this  information,  we  could  compute  a  magntitude  for  every  event,  this 
magnitude  being  def  ined  as  an  average  value  over  all  the  stations.  With  the  aim 
of  comparing  our  results  with  ISC  magnitudes  we  computed  the  mb(K)  -  mb(ISC)  (1)  and 
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also  mb(LDG)  -  mb(lSC)  (2). 

We  reported  these  differences  on  a  map  showing  the  event  locations  (see  Figure  4) . 

For  both  cases  (1)  and  (2)  the  differences  are  quite  similar  for  each  event.  One 
important  result  shown  on  the  map  is  the  geographic  partitioning  of  those  differences. 
We  noticed  that  for  the  events  located  in  the  middle  of  Shagan  River  test  site  all 
the  differences  were  zero  or  greater  than  0.1  while  the  events  located  in  the  North- 
East  part  of  the  test  site  show  negative  differences.  Below  we  will  relate  this 
feature  with  the  structure  at  the  test  site. 

In  conclusion  of  this  first  part  we  were  able  to  say  that  the  ’automatic  process'  we 
followed  to  obtain  a  magnitude  from  digital  data  provides  a  reliable  tool  to  compute 
magnitude  accounting  for  the  fact  that  the  'a'  value  should  be  used  and  is  character¬ 
istic  of  the  station. 

The  values  of  magnitudes  are  all  relative  since  we  did  not  have  any  reference  for 
the  yield. 


SPECTRAL  ANALYSIS 

We  were  interested  in  both  source-function  spectral  behavior  and  receiver  function. 

We  followed  a  method  described  by  different  authors  and  based  on  the  Mueller-Murphy 
(1971)  formula  for  the  theoretical  source-function. 

We  first  grouped  the  stations  into  six  subarrays  based  on  the  geological  features  of 
the  region.  For  each  subarray  the  region  was  assumed  geologically  homogeneous  (4 
stations  were  excluded  because  of  their  special  location  far  from  the  other  stations) . 
Then,  we  computed  the  wave  spectrum  from  the  spectral  density.  We  evaluated  the 
noise  spectral  density  from  a  window  10.24  s.  wide  ending  just  before  the  first  arrival 
and  the  signal  spectral  density  from  a  2.56  s.  wide  window  starting  at  the  onset 
of  the  first  arrival.  We  assumed  that  the  signal  was  the  convolution  of  a  source- 
function,  an  elastic  transfer  function  and  a  receiver  function. 

We  first  used  the  Mueller-Murphy  reduced  displacement  potential  for  an  underground 
nuclear  explosion  as  a  fource-function  to  deduce  the  transfer  function  for  the  path. 

We  based  this  calculation  on  absorption  band  model  described  by  Anderson  and  Given 
(1982).  We  computed  six  different  models  for  our  six  paths  from  the  source  area  to 
the  different  subarrays. 

With  these  results  and  the  data,  we  were  able  to  deduce  a  source-function,  in  theory. 
The  unknown  parameters  were  the  depth  and  the  yield,  but  the  inversion  program  did 
not  provide  a  good  results  because  the  parameter  variations  only  weakly  influence 
the  resulting  model. 

A  great  variation  in  the  yield  does  not  imply  a  great  variation  in  the  corner  frequency 
of  the  source-function.  Thus,  the  inversion  was  not  able  to  converge  to  a  well-defined 
and  unique  value. 
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TEMPORAL  ANALYSIS 

With  the  above  work,  we  can  define  the  limits  of  spectral  analysis  by  using  models. 

We  applied  the  same  method  but  in  the  time  domain  to  deduce  the  transfer  function 
by  using  the  Mueller-Murphy  model  as  a  reference  and  finally  a  source-function  for 
each  event. 

The  transfer  functions  computed  for  all  the  stations  defining  a  subarray  are  as  con¬ 
sistent  in  the  time  domain  as  in  the  frequency  domain,  except  for  the  subarrays 
located  in  the  Massif  Central  and  Morvan  region.  We  show  the  results  for  Vosges 
region  in  Figure  5. 

From  the  transfer  function,  we  deduced  receiver  function,  limited  in  that  we  only 
used  teleseisms  and  a  narrow  azimuthal  aperture. 

The  application  of  this  function  should  be  limited  to  thse  particular  events. 

Using  the  same  convolution  relation,  we  computed  source-functions  for  all  the  events 
from  the  recorded  signals  and  the  previously  computed  transfer  functions.  We  located 
all  the  source  function  on  a  map.  The  use  of  the  Mueller-Murphy  model  to  determine 
the  transfer  function  strongly  influenced  the  shape  of  the  resulting  receiver  function, 
but  we  could  notice  differences  in  the  shapes  of  these  source-functions.  A  site  effect 
was  indicated  in  the  North-East  part  of  Shagan  River  test  site  (see  Figure  6) .  There¬ 
fore,  we  divided  the  events  into  two  sets,  one  contained  the  events  located  in  the 
North-East  part  of  the  site  and  the  other  containing  the  events  located  on  the  middle 
part . 

We  computed  transfer  functions  again  but  using  only  the  second  set  of  data  to  avoid 
any  contamination  by  the  site  effect. 

The  results  were  not  very  different  in  either  temporal  response  or  frequency  content. 

We  noticed  only  that  the  introduction  of  the  first  set  of  data  in  the  computation  implied 
an  attenuation  in  the  spectral  shape  without  changing  the  content. 

New  source-functions  were  computed  from  these  "more  accurate"  transfer  functions.  That 
did  not  change  the  source-function  shape  for  the  events  in  the  central  part  of  the  test 
site,  but  a  second  peak  was  apparent  for  the  other  events  just  after  the  first  peak. 

In  order  to  confirm  the  validity  of  the  procedure  we  computed  the  "synthetic  signal" 
by  convolution  of  the  source-function  and  transfer  function  obtained  above. 

As  shown  in  Figure  7  for  two  stations,  we  found  a  very  good  agreement  between  this 
"synthetic  signal"  and  the  actual  signal. 

Then  we  tried  to  locate  the  anomaly  responsible  for  the  site  effect  using  a  representa¬ 
tion  of  the  ray  paths  for  the  P  waves  from  the  source  to  the  stations. 

We  first  plotted  the  ray  paths  near  the  stations  and  we  could  make  two  remarks: 

-  Taking  the  stations  with  extremum  azimuthal  angles  we  found  that  all  the  rays  were 
contained  in  a  beam  with  a  width  of  less  than  100  km  at  450  km  from  the  source  area 
which  means  that  the  paths  in  the  upper  mantle  are  about  the  same,  especially  in  depth 
regardless  of  the  station  and  source  locations. 

-  The  differentiation  of  the  signal  for  any  two  stations  was  already  clear  very  close 
to  the  source  area. 

Then,  we  plotted  the  ray  paths  around  the  test  site.  Since  the  events  were  very  close 
to  each  other,  we  could  see  some  overlappings  in  the  paths  from  two  different  events 
to  two  different  stations.  Therefore,  we  tried  to  compare  pairs  of  source-functions  in 
order  to  find  any  similitude  in  the  shape.  This  similitude  was  high  when  we  took  a 
pair  of  source-functions  which  overlapped  at  less  than  30  km  from  the  sources  as  shown 
on  Figure  8.  But  from  30  km  to  100  km  we  did  not  find  any  highly  correlated  pai-r  of 
source-functions  which  means  that  the  anomaly  was  very  close  to  the  source  area. 
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Because  the  events  did  not  cover  a  large  area  of  the  test  site,  we  were  not  able  to 
estimate  the  boundaries  for  this  anomaly  but  by  observing  the  evolution  in  the 
source-function  shape  in  Figure  9,  we  could  give  a  southern  limit  located  along  a 
Southeast-Northwest  trend. 

A  hypothesis  about  this  anomaly  is  that  the  rays  are  crossing  a  lens  with  a  very  slow 
velocity  generating  a  'secondary  source.'  This  hypothesis  seems  to  agree  with  our 
observations. 


CONCLUSION 

We  were  able  to  provide  a  reliable  method  to  compute  magnitude  from  digital  records 
with  some  assumptions  about  the  data.  This  method  was  at  least  as  accurate  as  the 
graphic  method. 

The  study  of  the  data  in  the  spectral  domain  allowed  us  to  determine  an  attenuation 
model  for  six  paths  from  the  source  area  to  the  subarrays.  Then  we  showed  that  these 
paths  could  be  characterized  by  a  transfer  function  applicable  to  the  event-France 
part  of  the  path  and  a  receiver  function  characteristic  of  each  station. 

From  these  transfer  functions,  we  computed  experimental  source-functions.  The  result 
showed  a  strong  site  effect  at  the  North-East  part  of  Shagan  River  test  site  which 
appears  in  the  source-function  as  a  secondary  arrival  probably  due  to  a  low  velocity 
region. 
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